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ABSTRACT 

Context. Mass accretion rate determinations are fundamental for an understanding of the evolution of pre-main sequence star 
circumstellar disks. 

Aims. Magnetospheric accretion models are used to derive values of the mass accretion rates in objects of very different 

properties, from brown dwarfs to intermediate-mass stars; we test the validity of these models in the brown dwarf regime, 

where the stellar mass and luminosity, as well as the mass accretion rate, are much lower than in T Tauri stars. 

Methods. We have measured simultaneously two infrared hydrogen lines, Pa;3 and Br.y, in a sample of 16 objects in the 

star-forming region p-Oph. The sample includes 7 very low mass objects and brown dwarfs and 9 T Tauri stars. 

Results. Brown dwarfs where both lines are detected have a ratio Pa/3/Br.yOf ~ 2. Larger values, ^3.5, are only found among 

the T Tauri stars. The low line ratios in brown dwarfs indicate that the lines cannot originate in the column of gas accreting 

from the disk onto the star along the magnetic field lines, and we suggest that they form instead in the shocked photosphere, 

heated to temperatures of ~ 3500 K. If so, in analogy to veiling estimates in T Tauri stars, the hydrogen infrared line fluxes 

may provide a reliable measure of the accretion rate in brown dwarfs. 



1. Introduction 

In pre-main sequence stars, circumstellar disks feed mat- 
ter onto their central stars over a period of a few million 
years. This accretion process does not significantly alter 
the properties of the star during most of this time; how- 
ever, it can have a significant influence on the lifetime 
and evolution of the disks and of the planetary systems 
that may form. The physical conditions of the accretion 
process, in particular the mass accretion rate, can be in- 
vestigated by studying the line and continuum emission 
produced by the accreting material. The model which has 
been more successful in explaining the accretion phenom- 
ena in T Tauri stars (TTS) is magnetospheric accretion, 
where the disk is truncated by the effect of a stellar mag- 
netic field near or within the corotation radius. The disk 
material, which is slowly drifting radially toward the star, 
reaches the truncation radius and is then lifted above the 
disk midplane and accretes onto the star along magnetic 
field lines, impacting on the star at approximately the 
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escape velocity. A shock forms on the stellar surface, cre- 
ating a hot spot which emits the excess UV continuum 
and lines observed in TTS. Optical and IR line emission 
is expected to come from the accreting columns of gas, 
where the temperature has to be of the order of 10000 
K. Magnetospheric accretion models have been developed 
in detail under the assumption of dipole magnetic field 
by, e.g., Hartmann et al. H1994fl . Calvet and GuUbring 
itTMSjl . MuzeroUe et al. l|1998lllj001|l . Lamzin ipg5..1998J , 
to predict the expected amount of veiling and line profiles 
and intensities. They have been used to derive accretion 
rates for TTS, brown dwarfs (BDs), some intermediate- 
mass TTS and Herbig Ae stars (Gullbring et al. [T^M 
MuzeroUe et al. IWl , lUU^, I^OMl VMM Natta et al.[ 
Calvet et al. ,;2Q0_4 Garcia Lopez et al. 



Testing magnetospheric accretion models is therefore 
very important and extensive checks have been carried out 
in several of the above papers. They have been mostly 
focussed on TTS, where it is possible to observe the 
accretion-driven emission over a large range of wave- 
lengths, both in lines and continuum. Although it is likely 
that the topology of the stellar magnetic field is much 
more complex than the dipole field for which most mod- 
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els have been developed, the models seem able to account 
satisfactorily for TTS activity (see Bouvier et al. H2(J(J6I 
and references therein). 

The situation is different for BDs, where at present the 
accreting matter is usually detected through optical and 
near-IR line emission. In fact, the vast majority of the 
existing estimates of Mace in BDs are derived by fitting 
the observed Hq, profiles with emission in magnetospheric 
accreting columns (Muzerolle et al. I2UU3L I2(JU5[ Natta et 
al. ED04j) or by means of secondary indicators, such as the 
luminosity of the hydrogen near-IR lines or of the Ca IR 
triplet, which have been calibrated using the Ha-measured 
accretion rates (Natta et al. |1T)04, Mohanty et al. 2005). 

This paper reports the result of a project aimed at 
testing the capability of current magnetospheric models 
to describe quantitatively the accretion in very low mass 
objects and BDs by comparing the observed ratios of two 
IR hydrogen recombination lines, namely Pa^ and Br-y, to 
the model predictions. These lines are seen in emission in 
most accreting TTS and BDs; their ratio can be measured 
relatively easily, and is not very sensitive to extinction. 
Muzerolle et al. (|2001|l models predict that for the rela- 
tivly large accretion rates of TTS (>10~* Mo/yr), both 
lines are optically thick and have a ratio Vap/Yii^'^ 4, 
typical of thermalized emission in a gas at about 10'' K. 
As Mace decreases below ^ 10"^ Mo/yr, first Br^ and 
then Pa^ become optically thin, and their ratio increases 
to much larger values. The sample of TTS studied by 
Muzerolle et al. f2001 ) does not have objects with very low 
Mace, and this particular prediction could not be tested. 
However, in our study of the accretion properties in p-Oph 
BDs (Natta et al. I2004II we found that two objects, with 
Maee'^W~^ Mo/yr, have 2, inconsistent with magne- 
tospheric accretion models. 

Unfortunately, the two lines had not been observed si- 
multaneously, and a chance low ratio due to time variabil- 
ity could not be ruled out. We present here the results of si- 
multaneous observations of Pa^ and Br-^ in a sample of 16 
objects in p-Oph, 7 of which are very low mass objects and 
BDs and the other 9 higher mass stars, added for compar- 
ison with the Taurus TTS. The observations and results 
are presented in Sec. 2 and 3, respectively. The reliability 
of the derived line ratios is discussed in Sec. 4. Comparison 
with the magnetospheric accretion models follow in Sec. 5. 
Sec. 6 contains summary and conclusions. 

2. Observations and data reduction 

2.1. Sample 

We have selected 16 objects in p-Oph which have been 
previously studied by Natta et al. IMiH 120061: 15 

of them are known to have Pa^ in emission, while in one 
case (ISO-Oph 033) Pa/3 was not detected, and we added 
it to our sample because it is the lowest mass object in 
p-Oph with a measured J-band spectrum. 

The adopted spectral type, effective temperature, lu- 
minosity and mass of the sample are given in Table 1. The 



spectral classification of objects in p-Oph is often very un- 
certain, as discussed, e.g., by Luhman and Rieke H1999|l . 
Wilking et al. p)05|l . Doppmann et al. H2005|l . Natta et 
al. (|200C3). For 7 very low mass objects and BDs (ISO- 
Oph 023, ISO-Oph 030, ISO-Oph 033, ISO-Oph 102, ISO- 
Oph 160, ISO-Oph 164 and ISO-Oph 193 ) we a dopt the 
stellar parameters derived by Natta et al. 1)2002(1 by com- 
paring low resolution simultaneous J, H and K spectra to 
template field stars and to model atmospheres. The un- 
certainties are of order ±150 K in Toff, ±0.2 dex in L*, 
±1 mag on Ay. 

For the other stars, we estimate the spectral type from 
the J-band medium resolution spectra obtained in this pa- 
per, as described in the Appendix. The corresponding Tcff 
is from Kenyon & Hartmann (1995) . The extinction A j is 
derived from the 2MASS (J-H) and (H-K) colors, by de- 
reddening the objects to the locus of the classical T Tauri 
stars defined by Meyer et al. 11997'), as discussed in Natta 
et al. H2006|l . We adopt an extinction law characterized 
by E(J-H)/E(H-K)=1.54 (Cardelh et aLUnHHl Kenyon et 
al. H1998|l . appropriate for Ophiuchus. The stellar luminos- 
ity is computed from the de-reddened J magnitude and the 
bolometric correction corresponding to the assigned spec- 
tral type (for a distance D = 150 pc). The uncertainties 
on Tcff and L^, are often large, ranging from ±100 K to 
±650 K in Tcff and between ±0.2 and ±0.5 dex in L^,. 

The location of the sample objects on the HR diagram 
is shown in Fig. 12. II together with the evolutionary tracks 
of D'Antona & Mazzitelli (|1997|l . used to estimate the 
object masses. The sample contains 7 objects with mass 
<0.1 Mq, and 9 with M^ in the range - 0.35 - 1 Mq. 

2.2. Near-infrared spectroscopy 

Near-infrared spectra in the J and K bands were obtained 
for all our targets using the ISAAC near-infrared camera 
and spectrograph at the ESO-VLT UTl telescope. The 
observations were carried out in Visitor Mode on May 19- 
20, 2005. We used the 0.6 arcsec slit and the medium res- 
olution grisni that offered a ~860 resolution across the 
wavelength 1.0-1.3 pm and ~750 at 1.95-2.45 pm. The 
observing sequence was composed of a set of pairs of spec- 
tra with the target in different positions along the slit, to 
allow for an efficient and accurate sky subtraction. The 
on-source times varied from 30 to 60 min depending on 
target brightness and observing band (the J band obser- 
vations were typically 30-40% shorter than the K band 
ones). Standard calibrations (flats and lamps) and tel- 
luric standards spectra were obtained for each observation. 
Wavelength calibration was performed using the lamp ob- 
servations. 

The Pap and Br-y lines were observed as "simultane- 
ously" as possible, to minimize the uncertainties due to 
the expected variability; to maintain a good temporal co- 
herence of the observations of the two emission lines, every 
target was observed in a time sequence Br-y-Pa/3-Br^. A 
complete observing cycle lasted about 2 hours. 
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Table 1. Observational and physical parameters of the p-Oph objects. Column 1: name of object; 2: spectral type; 3: 
J magnitude; 4: K magnitude; 5: J-band extinction; 6: luminosity; 7: mass; 8: effective temperature; 9: Pa^ equivalent 
width; 10: Br.^ equivalent width; 11: Pa^/Br-y flux ratio; 12: mass accretion rate. 



Obj. 


Spectral 


J 


K 




i* 


M. 


Teff 


Ew(Pa/3) 


Ew(Br7) 


R 


M 




Type 


(mag) 


(mag) 


(mag) 


(ie) 


(Mq) 


(K) 


(A) 


(A) 




(MQ/yr) 


ISO-Oph 002a 


MO 


12.838 


9.545 


3.1 


0.55 


0.4 


3850 


1.9±0.2 


2.1±0.3 


1.6±0.3 


-8.73 


ISO-Oph 023 


M7 


14.844 


12.143 


2.4 


0.04 


0.04 


2650 


1.2±0.2 


1.6±0.2 


1.7±0.4 


-9.69 


ISO-Oph 030 


M6 


12.57 


10.92 


0.9 


0.07 


0.06 


2700 


< 0.9 


< 6.0 


- 


<-9.54 


ISO-Oph 033 


M8.5 


16.45 


13.93 


2.2 


0.01 


0.015 


2400 


< 0.7 


< 2.7 




<-10.83 


ISO-Oph 037 


K5 


15.052 


10.224 


3.9 


0.16 


0.7 


4350 


4.4±0.5 


0.8±0.3 


3.8±1.5 


-9.61 


ISO-Oph 083 


K9 


12.257 


9.251 


3.0 


0.81 


0.4 


3900 


2.1±0.2 


< 0.9 


> 5.2 


-8.34 


ISO-Oph 102 


M6 


12.433 


10.766 


0.9 


0.08 


0.06 


2700 


1.4±0.3 


1.8±0.3 


2.0±0.5 


-9.17 


ISO-Oph 105 


K9 


12.547 


8.915 


4.0 


1.5 


0.35 


3900 


1.3±0.2 


< 3.0 


> 1.0 


-8.03 


ISO-Oph 115 


MO 


15.616 


11.486 


3.7 


0.07 


0.6 


3850 


1.1±0.4 


< 4.8 


> 0.3 


-10.81 


ISO-Oph 117 


K8 


13.325 


9.978 


3.1 


0.35 


0.6 


3950 


6.5±0.3 


3.3±0.1 


3.0±0.2 


-8.63 


ISO-Oph 155 


K3 


11.322 


7.806 


3.2 


3.1 


1.0 


4730 


1.3±0.1 


< 1.8 


> 0.9 


-8.28 


ISO-Oph 160 


M6 


14.148 


11.947 


1.9 


0.04 


0.05 


2700 


1.2±0.2 


2.3±0.4 


1.5±0.3 


-9.65 


ISO-Oph 163 


K5 


11.378 


8.271 


2.7 


1.6 


0.6 


4350 


3.3±0.2 


1.1±0.4 


4.6±1.7 


-7.89 


ISO-Oph 164 


M6 


13.27 


11.08 


1.9 


0.09 


0.06 


2700 


< 0.3 


< 2.8 




<-9.39 


ISO-Oph 166 


K5 


10.754 


8.464 


2.2 


1.9 


0.6 


4350 


1.7±0.2 


< 1.2 


> 3.9 


-8.19 


ISO-Oph 193 


M6 


13.611 


11.086 


2.3 


0.1 


0.06 


2650 


1.7±0.2 


1.9±0.4 


2.2±0.5 


-8.88 




3.8 3.6 3.4 3.2 

Log T,„ [K] 

Fig. 1. Location of the sample objects on the HR diagram. 
The evolutionary tracks (solid lines, labelled with the stel- 
lar mass) are from D'Antona and Mazzitelli H1997(l : dashed 
lines show isocrones at 0.5, 1 and 5 Myr, as marked. 

The spectra were reduced using standard procedures 
in IRAF. After checking that the two Br^ data sets were 
not significantly different (see §4.1), we combined all the 
exposures to obtain the final spectrum. The portion of the 
resulting spectra centered around the Pa^ and Br-^ lines 
is shown in Fig. 12.21 



3. Results 

The equivalent widths of Pa^ and Br.^ were determined by 
fitting gaussian profiles to the spectra using standard iraf 
routines. The uncertainties were estimated from the uncer- 
tainties on the continuum location and from the equivalent 
width of noise spikes. We have checked that the equivalent 
width obtained by integrating the observed flux in excess 
of the continuum in a fixed wavelength range does not 
deviate from the values of the gaussian fit by more than 
la. 

In 13 objects Pa^ is clearly seen in emission, with 
equivalent width between ~ 1 and 6.5 A. In 3 cases, the 
line is not seen and we report in Table 1 ia upper lim- 
its. In all objects, the emission lines are not resolved, with 
the exception of two of the sample TTS (ISO-Oph 163 and 
ISO-Oph 155), where we clearly detect red-shifted absorp- 
tion. However, the resolution of our spectra is not suited 
for a study of the line profiles. Br..^ is generally weaker, 
and we have detected it in 8 objects only; in some cases, 
the equivalent widths have large uncertainties. 

A comparison of the equivalent width of this paper 
with the measurements reported by Natta et al. H2004L 
12006(1 shows the importance of measuring line ratios from 
spectra obtained simultaneously, as variations of factor 
of 2 are possible. Of the three objects with non de- 
tected Pa/3, in two cases the current upper limits are con- 
sistent with our previous observations. In one case (ISO- 
Oph 164), however, Pa^g was clearly detected by Natta et 
al. 1)2004(1 . with an equivalent width of 0.8A, significantly 
larger than our current upper limit (0.3 A). 

We compute line fiuxes from the equivalent widths and 
the flux of the nearby continuum. We do not correct the 
equivalent width for underlying photospheric absorption, 
since the expected equivalent widths are small (<0.5 A, 
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Fig. 2. Pa/3 (left panel) and Br7 (right panel) spectra 
of the target objects. The observed fluxes have been nor- 
malized to the continuum near the line and shifted for an 
easier display. 



Wallace ate al. HUU^ for objects with Toff <5000 K. As 
it was not possible to flux-calibrate the spectra, we esti- 
mated the continuum flux using the 2MASS values of the 
J and K magnitudes, proceeding as follows. The spectra 
were integrated over the 2MASS filter profiles and scaled 
to the measured 2MASS values. The resulting line fluxes 
were then corrected for extinction, according to the ex- 
tinction law of Cardelli et al. I1H8HI for i?y=4.2 and the 
values of Aj given in Table 1. 

The ratio of the line fluxes Pa^j/Br-y {R in the follow- 
ing) is given in Table 1, which gives also the la errors on 
i?, computed taking into account only the errors on the 
measured equivalent widths and on the 2MASS photom- 
etry. They range from 3 to ~ 40%, in objects with very 
weak Br J. The possible systematic errors introduced by 
our flux calibration method and the extinction correction 
are discussed below. 

Table 1, last column, gives the mass accretion rate 
Mace, computed from the luminosity of Pa,p, using the 



relation between L(Pa^) and the accretion luminosity de- 
rived by MuzeroUe et al. (|1998|l and Natta et al. (fM12| : 

log Lace/ Lq ^ IM log L{Pp)/Lq+ 4 (1) 
The mass accretion rate is then computed from hacc- 

Mace ^ Lace R*/iGM,) (2) 

The uncertainties on Mace are dominated in most cases 
by uncertainties on the stellar parameters, and by the scat- 
ter of eq.(l) for very low luminosity objects. They can 
easily be of the order of one dex in Log Mace- 

4. Discussion 
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Fig. 3. R is shown as a function of (Panel a) and of 
Mj, (Panel b) for the p-Oph objects; dots are measured 
values, arrows lower limits. Crosses plots R values for the 
Taurus TTS from MuzeroUe et al. ifDMjl . 

The values of R span a rather large interval, varying 
from less than 2 to >5. We plot the values of R in Fig.Oas 
a function of the luminosity of the central object and of its 
mass. Inspection of Fig. 13 shows that the 4 low luminos- 
ity (L*< O.ILq), low mass (M*< O.O6M0) objects with 
measured R (ISO-Oph 023, ISO-Oph 102, ISO-Oph 160, 
ISO-Oph 193) have R<2. Values of R>3.5 are only found 
in more massive objects, with M*> 0.4 M0. 

In flg. 131 we have added to our objects the TTS 
in Taurus for which the ratio R has been measured by 
MuzeroUe et al. H1998|l . For better consistency, we have re- 
computed the R values from the line equivalent widths and 
A J given by MuzeroUe et al. H1998|l . and 2MASS photome- 
try. The differences with the values reported by MuzeroUe 
et al. (|1998|l are generally small. 

The Taurus sample does not include any very low mass 
object and is limited to stars with M*>0.25 Mq, but it 
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complements nicely our p-Oph data for TTS. In particular, 
one can see how, on one side, large values of i?(>3.5) are 
observed only in higher mass objects, while on the other 
all the observed BDs have R<2. 

There are few higher mass objects with 2. One 
is ISO-Oph 002, which will be discussed in Sec. 5. For the 
three Taurus objects with similarly low R (GI Tau, FS Tau 
and DG Tau), we note that the Taurus observations of Pa^ 
and Br^ are not simultaneous, and line variability may af- 
fect by large factors the R value of individual objects. For 
example, the Pa^ equivalent width measured for GI Tau at 
two different epochs by Folha and Emerson 1)2001(1 varies 
by more than a factor of 2. 

The low values of R are particularly interesting, as we 
will dicuss in the following. However, before doing that, 
we need to examine the potential sources of error in the 
estimates of R, which may undermine the significance of 
our results. 



4.1. Time variability 

Line and continuum variability are well known in pre-main 
sequence stars of all masses. We have measured the two 
lines simultaneously, using an observational scheme Br^- 
Pa/3-Br-y intended to minimize the effect of variability on 
R. In fact, it turns out that on the time scale of 1-2 hours, 
spanned by our observations, variability is not large, and 
the Br^ equivalent widths measured before and after Pa^ 
are the same within the errors. 

Note, however, that it has not been possible to cal- 
ibrate our spectra photometrically, so that, in order to 
compute the line ratio i?, we had to rely on magnitudes 
from the literature. 2MASS observations have been per- 
formed simultaneously in the three bands J,H,K, so that 
the values of R assigned to each object can be consid- 
ered a correct, albeit "snapshot", value, as long as there 
are no strong variations in the (J-K) colors of the stars. 
This seems to be indeed the shown by Fig. 0] 

which compares the values of R obtained by using the J 
and K magnitudes measured by DENIS and by Barsony 
et al. p997|l . respectively, to those in Table 1, computed 
with 2MASS photometry. All sets of photometric obser- 
vations have been obtained simultaneously in the J and 
K bands. The differences are always within 2a\ the five 
values of R lower than ~ 2 remain low, while the R values 
of the order of 4 or larger remain so. 



4.2. The extinction correction 

The values of R depend on the differential extinction be- 
tween the J and K band. This extinction correction can 
be large in p-Oph objects, which are often very embedded, 
and introduces an additional uncertainty on R. In partic- 
ular, underestimating Aj will result in values of R lower 
than the true ones, since R is proportional to ' (for 

the adopted extinction law). Therefore, to push observed 



E- 

CD 



Q 



CO 



I*..- 




R (Pa/Br^) [2MASS] 



Fig. 4. Comparison of R using three different sets of near- 
IR photometry to derive the continuum flux at the line 
wavelengths. The horizontal axis shows R computed from 
the 2MASS photometry; on the vertical axis, the results 
obtained using DENIS (squares) and Barsony et al. H1997|) 
observations (triangles). 

values of R from ~ 2 to ^ 4, one needs to increase the J 
band extinction by 1.37 mag. 

Such an increase of Aj with respect to the values given 
in Table 1 is not consistent with the J, H, K spectra 
and photometry of the 4 low mass, low R objects (ISO- 
Oph 023, ISO-Oph 102, ISO-O ph 160 and ISO-Oph 193) 
studied by Natta et al. I)2002|l . who estimate uncertain- 
ties in A J of ^^0.3 mag. In fact, independent estimates of 
Ay for two objects (ISO-Oph 030 and ISO-Oph 102) by 
Wilking et al. H2()()5|l result in values of Aj even lower 
than those adopted in this paper. 

It is more difficult to estimate the Aj uncertainties for 
the 8 stars in our sample for which Aj is derived from 
the (J-II)-(II-K) colors. Luhman and Rieke (jl999|l . using 
the Barsony et al. (1997) photometry and a more complex 
scheme than the simpler one we have implemented, derive 
values of Aj for the 5 objects we have in common within 
0.3 mag of our estimates. In one case (ISO-Oph 166) there 
is also in the literature an estimate of the extinction from 
optical spectra (Wilking et al. I2005|l . which agrees with 
our value within 0.1 mag. 

In summary, although it is difficult to rule out the pos- 
sibility of a large error on Aj in some individual object, 
we think it is unlikely that this can systematically affect 
our results. In particular, it does not seem possible to in- 
terpret all the low R values as due to an underestimate of 
the extinction. 

5. Comparison with magnetospheric accretion 
models 

The main aim of this work was to extend toward very 
low mass objects the comparison of the observed ratios 
Pa/j/Br^ to the prediction of magnetospheric accretion 
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models. In these models, gas is channeled from the disk 
onto the star along the lines of the stellar magnetic field 
and a shock forms where the accreting matter impacts on 
the stellar photosphere. The accretion luminosity is emit- 
ted mostly by the heated photosphere below the shock (the 
hot spots on the stellar surface), with the rest coming out 
from the pre-shock region. In accreting TTS, the heated 
photosphere and the pre-shock region are responsible for 
the continuum which veils the optical spectrum and for 
the strong ultraviolet continuum and line spectrum seen 
in these stars. The optical lines are believed to form not 
in the shocked region, which is very optically thick, but in 
the accreting columns, heated to roughly uniform temper- 
atures (of order 8000 - 12000 K) not by the shock emission 
but by some alternative mechanism, possibly wave heat- 
ing. These models have been very successful in explaining 
the complex phenomenon of TTS activity (e.g., Bouvier 
et al. I200(il and references therein) . 

In the magnetospheric accretion models, for fixed stel- 
lar parameters and field geometry, R depends on the mass 
accretion rate Mace and temperature of the accreting gas. 
As mentioned, the gas temperature is not computed self- 
consistently, but has been constrained by fitting different 
lines of various atomic species (see Muzerolle et al. l2()()l|l . 
Within these constraints, the infrared lines are thermal- 
ized for the gas densities expected for average values of 
Mace- Both the line luminosities and R do not change 
with decreasing density until A/qcc^SIO^^ Mo/yr, at which 
point the line optical depth goes below unity first at Br^ 
and then Pa^ and R quickly increases to very large values. 
There are additional free parameters, namely the exact 
location of the base of the accreting region and its incli- 
nation to the line of sight, which, however, do not change 
this basic picture. 

The predictions of magnetospheric accretion models 
are compared to our results in Fig. El which shows that 
the models reproduce quite well the bulk of the R values 
in the Taurus TTS (Muzerolle et al. l^niIT|l . 

The situation is different for the p-Oph BDs, which 
have much lower Mace than TTS (e.g., Muzerolle et 
al. [2003, 2003 Natta et al. IMHl Mohanty et al. IMISll . 
Fig. shows model prediction for stellar parameters typi- 
cal of BDs and Macc^ 10~^ Mo/yr. The gas temperature 
required to fit intensity and profiles of Hq is T>8000 K. 
In these models R ^ 3.5 for Mace— 10~^ MQ/yr, and in- 
creases rapidly with decreasing Mace- It is clear that they 
cannot reproduce the observed low R values of the p-Oph 
BDs. 

2 is obtained if the lines form in a cold, optically 
thick gas. This is shown in Fig. El which plots the R val- 
ues expected as a function of the gas temperature T in 
different physical situations. The line at i?~ 6 shows the 
predictions of Case B, i.e., the value of R when the lines 
are optically thin and the level population is dominated 
by radiative cascade from the continuum. If the hydrogen 
levels from which the near-IR lines originate {n > 3) are 
in LTE, the dependence of R on the temperature is differ- 



ent for optically thin and optically thick lines. When the 
lines are optically thin, i?>8 for temperature T ~ lOOOOK 
and increases for decreasing T (dashed line). On the con- 
trary, when the lines are optically thick, R decreases for 
decreasing temperature, from ^ 4 for T> 8000 K to ^ 2 
for T - 3500 K. 

If indeed the lines are thermalized, the size of the emit- 
ting region can be estimated from Fig. which plots the 
luminosity of Pa/3 as a function of R for the Ophiuchus 
and Taurus sample. The three curves show the location of 
thermalized, optically thick IR line emission for different 
values of the projected area of the emitting region, from 
2 X 10^2 to 2 X 10^0 cm2 and a FWHM fine width of 100 
km s-i. The four BDs (ISO-Oph 23, ISO-Oph 102, ISO- 
Oph 160 and ISO-Oph 193) with i?~ 2, are well fitted by 
black body emission at T^ 3500 K and area 2 x 10^^ 
cm^. This is a small fraction of the projected photospheric 
surface, varying between 6 and 15% in the four objects. 
Note that at the resolution of our data (~ 350 km/s) the 
BDs lines are not resolved; as the V&p flux is proportional 
to the product of the emitting area times the line width, 
larger or smaller values of the latter will result in a propor- 
tionally smaller or larger values of the emitting surface. 

These properties (i.e., temperature higher than Teff, 
small covering fraction, high optical depth) remind those 
of the hot spots caused by magnetospheric accretion on 
the stellar surface of classical TTS. 

Models of the accretion shock have been computed by 
Calvet and Gullbring H1998|l for the parameters typical 
of TTS. They showed that the emission originates from 
three different regions, the pre-shock accreting columns, 
the post-shock gas and the heated photosphere below the 
shock (the hot spot) . The physical properties of these re- 
gions are controlled by the energy flux of the accretion 
flow JT: 



1- no 1 nlO -2 -1 / Mace N 

^ = 9.8x10 ergcm s ( ,o_SM^yr-0 

(l^y (3) 
VO.5M0/ V2R0/ Vo.Ol/ ^ ^ 

where / is the fraction of the stellar surface covered by 
the accretion spots. In TTS J- ranges between 3 x 10^^ 
and 10^2 erg cm"^ s~^; with these values of J-, the heated 
photosphere is very optically thick, and emits a continuum 
spectrum with temperatures much higher than Tea, which 
accounts for the UV veiling. The models are quite com- 
plex and cannot be simply scaled to the BD parameters. 
However, we note that the energy flux of the accretion 
flow in the four p-Oph BDs is ^ ~ 5 x 10^ — 5 x 10^ erg 
cm^2 s^^, much lower than in TTS (see Fig. From 
the Calvet and Gullbring models, one expects that both 
the enhancement of the temperature with respect to the 
undisturbed photosphere, and the depth of the affected 
region decreases with JT, so that it seems possible to ob- 
tain line emission with the required properties from the 
heated photosphere. 
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Fig. |S1 plots the observed values of i? as a function of 
T for Taurus and p-Oph objects. For Taurus, we have 
computed J- from the accretion rates, stellar parameters 
and the surface coverage of the accretion columns given 
by Calvet and Gullbring H1998|l . using their eq.(ll). For 
the p-Oph objects, we take the temperature and the area 
of the emitting region from the location of the object in 
Fig 13 assuming black body emission. This is possible for 
the 6 objects with i?<3.5. One can see that 5 of them (the 
4 BDs and ISO-Oph 002), which have R-^ 2, also have 
J- much smaller than the Taurus TTS. The sixth object 
(ISO-Oph 117) has i?~ 3.5 and T similar to the lowest 
values of Gullbring & Calvet. The correlation between R 
and J- stops at ^ 10^"; for higher T , all the objects 
have similar values of R. 

We propose a scenario where for high Mace the opti- 
cal and infrared continuum emission of the heated photo- 
sphere is optically thick and all the hydrogen lines (optical 
and IR) observed in emission originate from the accreting 
columns of gas. When Mace decreases below a threshold 
value the continuum emission of the hot spots becomes op- 
tically thin and the lines appear in emission; at the same 
time, the density of the accreting gas decreases and its line 
spectrum gets weaker. As a consequence, lines of lower 
optical depth have an increasing contribution from the 
shocked regions. In BDs, this region dominates the emis- 
sion of the near-IR hydrogen lines, while optical lines, such 
as Hq, still originates mostly in the accreting columns of 
gas. If this will be confirmed by further observations and 
by model calculations, the IR hydrogen lines may provide 
a measurement of Mace in BDs in analogy to continuum 
veiling in TTS, less model-dependent than the Hq profiles 
used so far. 

A similar picture is suggested by Edwards et al. (|2[)()()(l , 
who obtained high resolution profiles of Fa-,, in classical 
Taurus TTS. They find that the Pa^ width correlates with 
the veiling at 1 /im, which may be a tracer of Mace- The 
width of lines forming in the accreting gas columns should 
not depend on Mace for fixed stellar and disk parameters, 
and Edwards et al. argue that their results can be under- 
stood if the relative contribution of the accretion shock to 
the Pa^ emission increases as Mace decreases. 

6. Summary and conclusions 

This paper presents the results of simultaneous measure- 
ments of Pa^ and Br-,, in a sample of 16 objects in the 
p-Oph region. We measure the ratio of the line fluxes 
i?=Pa^/Br.y in 8 objects, 4 of which are spectrally con- 
firmed very low mass stars or brown dwarfs (BDs in this 
paper). In 5 other cases we could not detect Br-y, and we 
estimate lower limits to R. 

The 4 BDs all have i?~ 2; of the other stars, one has 
R^ 2, while the others have R> 3.5, similar to the val- 
ues measured in accreting TTS in Taurus (MuzeroUe et 
al. I10f)8|l . We discuss the measured R values in the con- 
text of magnetospheric accretion models where hydrogen 
recombination lines form in the accreting columns of gas 
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Fig. 5. iZ as a function of Mace for the p-Oph objects; 
dots show actual measurements, arrow lower limits. Mace 
has been derived from the Pa^ luminosity, as described 
in the text. Crosses show the location of Taurus TTS 
(from Muzzerolle et al. I1998|l . The different lines show 
the predictions of magnetospheric accretion models which 
roughly fit the UV veiling and the H„ profiles. Solid and 
short-dashed lines show models for TTS (M^=0.5 M©, 
R*=2 Rq) with inner and outer magnetosphere bound- 
aries at the disk of 2.2-3 R* (solid) and 2.8-3 R* (short- 
dashed). The temperature of the accreting gas increase 
from 7000 K (for Mace^ 10"^ Mg/yr) to 12000 K for 
Mace=\^^'^ M0/yr(see Muzerolle et al. lMITI for details). 
The other two lines show models for typical BD parame- 
ters (M^,=0.05 M0, R^=0.5 R©), magnetosphere bound- 
aries 2.2-3 R^, temperature of 12000 K (long-dashed) and 
10000 K (dot-dashed), respectively. 



that connect the disk to the star. The gas is heated to 
temperatures of ^ 8000 — 12000K, and the predicted val- 
ues of R are in all cases > 3.5. The low R observed in the 
p-Oph BDs cannot be explained by these models. 

We propose that, at the low Mace typical of the BDs, 
the hydrogen near-IR lines form instead in the shock- 
heated stellar photosphere, and discuss how it is possible 
to produce in these regions the required low temperature 
and high optical depths in the line (but not in the con- 
tinuum) that will result Pa^g and Br^ emission with in 
i?~ 2. This possibility needs to be investigated with the 
help of models similar to those of Calvet and Gullbring 
l)1998|l . which compute the emission of the various com- 
ponents of magnetospheric accretion (accreting column, 
shock-heated photosphere and post-shock gas) for low 
Mace values and the parameters typical of BDs. 

Finally, we note that the luminosity of the hydrogen IR 
recombination lines has been extensively used to derive the 
accretion rate in large sample of objects, especially when 
UV veiling could not be measured (e.g., Natta et al. 12006(1 . 
The calibration of the line luminosity - accretion relation- 
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Fig. 6. R is plotted as a function of the gas temperature 
Tgas for three textbook cases. The sohd hne shows the 
case of optically thick lines with LTE level population; 
the dashed line the case of LTE, optically thin lines. The 
dash-dotted line the value of R for case B recombination. 



ship has been derived empirically, using objects with Lqcc 
measured independently from veiling or (for very low mass 
objects) from the Hq, profiles, not from models. The results 
of this paper, which suggests that the IR lines may not 
necessarily form in the accreting columns, do not weaken 
their reliability as a tool to measure the accretion rate. In 
fact, if they form in the accretion spots, they may provide 
a measurement of Mace in BDs in analogy to continuum 
veiling in TTS, less model-dependent than the Hq, profiles 
used so far. 

Appendix A: Spectral classification of the 
Ophiuchus objects 

Many of our targets do not have a reliable spectral classi- 
fication in the literature, others have been classified using 
a variety of different methods, from low resolution opti- 
cal, to K-band or high resolution infrared spectroscopy. 
We attempted to use our J-band spectra to estimate in 
a uniform fashion the spectral types of all objects in our 
sample. The J-band was chosen because at this wavelength 
the photosphere dominates the emission even for Class II 
objects and the effect of veiling from the disk infrared ex- 
cess emission is thus limited, especially in the BDs. 

We computed the equivalent widths for a number of 
photospheric lines, following Wallace H2()()()|l . Due to the 
lower resolution of our spectra, we needed to adjust the 
limits for the computation of the equivalent widths given 
in Table 3 of Wallace et al.; thus, we recomputed the val- 
ues of the equivalent widths for all spectra in their list. We 
also computed the equivalent width values for all p-Oph 
Class II objects in our Natta et al. I)20(J6|I sample with a 
solid spectral classification in the literature, mainly from 
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R (Pa/Br^) 

Fig. 7. Pa^ flux versus R for Ophiuchus (dots and ar- 
rows) and Taurus (crosses). The three curves ending at 
i?'^ 4 show the LTE, optically thick line emission case. 
Along each line, the temperature increases from 1000 K 
to 15000 K, in steps of 1000 K (tick marks); the first mark 
corresponds to T=2000 K on the top curve, to T=3000 
K for the middle curve and to T=5000 K for the low- 
est one. Note that the values of R depend on T only. Each 
curve corresponds to a different value of the emitting area, 
2 X 10^^ cm^ (solid hne), 2 x 10^^ cm^ (long-dashed hne) 
and 2 x lO'^'^ cm^ (dotted line) from top to bottom, re- 
spectively. In all cases, we have assumed a Pa^j FWHM of 
100 km s~^; the Pap flux is proportional to the emitting 
area times the line width. Magnetospheric accretion mod- 
els (same as in Fig.El are shown by the roughly horizontal 
lines (same as Fig.lSJ. Slightly lower values of the gas tem- 
perature in the accreting columns will give lower value of 
F(Pap), while still roughly fitting the observations. Note, 
however, that in no case R< 3. 



Luhman & Rieke lTM^ or from Natta et aL l^UH^ To im- 
prove the possibility of a good spectral classification at late 
types we also defined a spectral index based on the water 
absorption features at the red edge of the band, in a simi- 
lar fashion as Testi et al. H2001fl . This index was computed 
only for the p-Oph Class II objects with known spectral 
type. In Figure lA.ll we show the run of the equivalent 
widths of [AlII] and [KI] and the water index as a func- 
tion of the known spectral type. For late spectral types 
(K-M) a combination of the various indices may allow a 
reasonable classification within a few subclasses, at ear- 
lier spectral types the correlations of the indices with the 
spectral type flattens significantly and the classification is 
more uncertain (if at all possible) . 

Using the equivalent widths and the water index, for 
each of the stars in our sample we derived a range of pos- 
sible spectral types compatible with the measured values. 
These are reported in Table lA.ll along with the spectral 
types found in the literature. In general our estimate of 
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Fig. 8. R values plotted against the energy flux of the 
accretion flow J^. Estimates for the p-Oph objects (dots) 
have been derived from their measured Mace, stellar pa- 
rameters and an emitting area estimated by tring heir lo- 
cation on the F(Pa/3)-i? with the predictions of LTE, op- 
tically thick emission (Fig.|7I); two objects (ISO-Oph 037 
and ISO-Oph 163) with i?>4.5 no estimate of T was pos- 
sible. The values for Taurus TTS have been computed 
from Mace and stellar parameters given by Muzerolle et 
al. H1998|l and the fraction fo stellar surface covered by the 
accretion spots computed by Calvet and GuUbring (1998|. 
In all cases we have used eq.(ll) of Calvet and Gullbring 



the spectral type is consistent with the values found in 
the literature. Following this comparison, for each target, 
we adopt a fiducial spectral type that is used through- 
out this paper to derive the stellar parameters (Column 
9). For the very low mass stars and brown dwarfs stud- 
ied in Natta et al. H2002|l we adopt their spectral type as 
our spectral type range is essentially based on the water 
index and is more uncertain than the classification made 
using the complete near infrared spectrum. For these ob- 
jects the adopted type is always within the range we derive 
here, with the exception of pOph-ISO 193, for which an 
earlier type (by 2 subclasses) is possible. For the stellar 
objects, the adopted spectral type is at the center of the 
range derived here, and it is usually consistent with the 
estimates in the literature. Some of the objects, however, 
have very uncertain spectral classification, this is reflected 
in the larger uncertainties in the effective temperature es- 
timate used in Sect. 12. II and Fig. 12.11 
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Fig. A.l. Run of the main spectral indices used for the 
spectral classification as a function of the spectral type. 
The top panel shows an index that measures the water 
vapour absorption at the red end of the J-band, the middle 
and bottom panels show the equivalent width of the [Kl] 
and [AlII] lines within the J-band. Open squares, triangles 
and diamonds show the values of the [KI] and [AUl] equiv- 
alent widths computed for the giants (luminosity class III), 
dwarfs (V) and subgiants (IV) of WaUace et al. (2000), re- 
spectively. Crosses show the values of the same indices for 
the objects with known spectral type in Ophiuchus. The 
solid lines are smoothed averages of the measured points, 
the dotted lines show the dispersion within the smoothing 
interval (±3 subclasses); they have been used to estimate 
the spectral type range for our target stars from the mea- 
sured indices. 



Table A.l. Line indexes and derived spectral types. Column 1: name of object; 2, 3, 4: H2O, KI and All indexes measured from this paper spectra; 5, 6, 7: 
same from Natta et al. ()2()()f)|l J-band spectra; column 8: spectral type range; column 9: adopted spectral type; column 10-13: spectral types from Natta et 
al. (|2(JU2|I . Luhman and Rieke H1999|l . Wilking et al. I)2U(J5|I and Doppmann et al. ()2005|l . respectively; column 14: other source names. 



Name 


H2O 


K I 


Al I 


H2O 


K I 


Al I 


ST 


ST 


ST 


ST 


ST 


ST 


Other Names 




This Paper 




Nea06 




range 


adopted 


Nea02 


LR99 


Wea05 


Dea05 




ISO-Oph 002a 


0.95 


1.9 


2.1 


0.88 


1.8 


1.6 


K8-M2 


MO 


- 


- 


- 


- 


B 162538-242238 


ISO-Oph 023 


0.75 


1.5 


0.5 


0.78 


2.7 


1.1 


M4-M9 


M7 


M7 


— 


— 


— 


SKSl 


ISO-Oph 030 


0.79 


2.2 


1.2 


0.82 


2.2 


1.4 


M4-M8 


M6 


M6 


M5-M6 


M5.5 


_ 


GY5 


ISO-Oph 033 


0.71 


2.7 


0.7 


0.73 


3.4 


0.5 


M7-M9 


M8.5 


M8.5 


_ 


_ 


_ 


GYll 


ISO-Oph 037 


0.99 


1.0 


1.9 


0.92 


1.7 


2.5 


K0-K9 


K5 


- 


K-M 


- 


K9 


LFAM3/GY21 


ISO-Oph 083 


0.93 


1.9 


2.0 


0.92 


1.7 


1.4 


K8-K9 


K9 


— 


— 


— 


— 


BL162656- 241353 


ISO-Oph 102 


0.81 


2.1 


1.1 


0.80 


2.7 


1.1 


M3-M6 


M6 


M6 


— 


M5.5 


— 


GY204 


ISO-Oph 105 


0.93 


1.5 


1.6 


0.92 


1.1 


1.8 


K8-M0 


K9 


- 


K5-M2 


- 


M3 


WL17/GY205 


ISO-Oph 115 


0.85 


0.8 


1.2 


0.88 


1.0 


2.1 


K9-M2 


MO 










WL11/GY229 


ISO-Oph 117 


0.93 


1.4 


1.3 


0.86 


1.2 


2.4 


K2-M3 


K8 










GY235 


ISO-Oph 155 


0.96 


0.9 


1.4 


0.94 


0.8 


1.3 


K0-K5 


K3 




K5-M2 






GY292 


ISO-Oph 160 


0.81 


2.2 


1.2 


0.81 


2.6 


1.4 


M3-M8 


M6 


M6 








B162737-241756 


ISO-Oph 163 


0.96 


1.2 


1.6 


0.91 


1.3 


1.8 


K0-K8 


K5 




K5-M2 






IRS49/GY308 


ISO-Oph 164 


0.70 


3.0 


1.0 


0.71 


3.5 


0.7 


M7-M9 


M6 


M6 




M3-M5 




GY310 


ISO-Oph 166 


0.93 


1.6 


2.1 


0.91 


1.4 


1.8 


K2-K8 


K5 




K6-M2 


KO 




GY314 


ISO-Oph 193 


0.86 


1.4 


1.3 


0.87 


1.4 


0.8 


M1-M4 


M6 


M6 








B162812-241138 
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